Adiabatic Phase Diagram on Degenerate Fermi 
Gas with Feshbach-Resonance 



S. Watabe, T. Nikuni*, N. Nygaardt, J. E. Williams*, and C. W. Clark* 

Department of Physics, Graduate School of Science, University of Tokyo, 

Tokyo, Japan 

* Department of Physics, Faculty of Science, Tokyo University of Science, 

Tokyo, Japan 

^Danish National Research Foundation Center for Quantum Optics, 
Department of Physics and Astronomy, University of Aarhus, 
DK-800 Arhus C, Denmark 
* Electron and Optical Physics Division, 
National Institute of Standards and Technology, Gaithersburg, 
MD 20899-8410, USA 



We determine the adiabatic phase diagrams for a resonantly- coupled system 
of Fermi atoms and Bose molecules confined in a harmonic trap by using 
the local density approximation. The key idea of our work is conservation 
of entropy through the adiabatic process. We also calculate the molecular 
conversion efficiency as a function of the initial temperature. Our work helps 
to understand recent experiments on the BCS-BEC crossover, in terms of the 
initial temperature measured before a sweep of the magnetic field. 

1. Introduction 

Degenerate Fermi gases with controllable interaction have been realized 
in resent experiments by making use of a Feshbach resonance. By ramping a 
magnetic field across the resonance, one can observe crossover between BCS 
and BEC superfluidity in a Fermi gas P^l Williams et al. 3 calculated equilib- 
rium phase diagrams for an ideal gas mixture of fermionic atoms and bosonic 
molecules, which provide qualitative understanding of experimental data of 
Refs Jll2l In this paper, we extend the work by Williams et al. 3 to include 
the resonant interaction using mean-field theory Calculating the entropy as 
a function of the temperature and the resonance energy, we show the paths 
of constant entropy traversed in the conventional phase diagrams as the res- 
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onance energy is varied adiabatically. On the basis of the adiabatic path, we 
determine the adiabatic phase diagrams against the resonance energy and 
entropy, labeled by the initial temperature of the Fermi gas measured before 
a sweep of the magnetic field. We also calculate the production efficiency 
of molecules during the adiabatic sweep of the magnetic field, by extend- 
ing the calculation given in the recent work by Williams et al 4 to include 
mean-field interactions. Our study gives an intuitive understanding of recent 
experiments on the BCS-BEC crossover using the adiabatic sweep process. 



2. Equilibrium Theory and Phase Diagrams 

The starting point of our theor y is th e coupled boson-fermion Hamilto- 
nian with two-component Fermi gasP^ Our mean-field calculation closely 
follows the procedure described in detail in RefEl The resulting grand canon- 
ical potential Q is given by 

O = Jdv ([e res - 2/.(r)]|$(r)| 2 + I|^L_l n {l - -e ves -2^ r ) ]} 
+ /(^{ [£ k ) -^ r )-^(r)]-|ln[l + e-^)]}). (1) 

The chemical potential \i is introduced to impose the constraint that the 
number of particles be conserved. The effect of a harmonic trap is included 
by the local density approximation, 8 which adopts the effect of local potential 
by replacing the chemical potential /i with /i(r) = fi — V a (r), where V a (r) 
is the harmonic trap for fermionic atoms, is the kinetic energy of a 

molecule, and is the kinetic energy of an atom. The local quasiparticle 

excitation energy is i?k( r ) = y [^k M r )] 2 + I A(r)| 2 , where the local gap 
energy is |A(r)| = tt|$(r)|, <3>(r) is the condensed molecular wavefunction, 
and k is the coupling constant of the resonance interaction. Finally, £res is 
the bare resonance energy, which is controlled by the magnetic field. 

We regard the equation of the number of particles and the gap equation 
as the simultaneous equations for the chemical potential fi and the local 
gap |A(r)|, for a given temperature T and a resonance energy e res . In this 
paper, we consider the case of narrow resonance. We set the bare coupling 
strength to a = yfpK — 0.4£p, where p is the peak density and sy is the 
Fermi energy of a pure gas of fermionic atoms. To avoid the divergence in 
the gap equation, we impose the gaussian cut-off. 

The condensed pair number N c is composed of the number of condensed 
molecules N mc and Cooper pairs N p . In FigQ] (A), we plot the condensed 
molecular fraction rj mc — N mc /N tot against the temperature T and the 
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Fig. 1. Phase diagrams of the condensed molecular fraction rj mc (A) and the 
Cooper pair fraction rj p (B). 



resonance energy e res . We also plot the condensed Cooper pair fraction 
rjp — Np/Ntot, which was defined by Yang^, in FigE](B). The energy scale is 
normalized by the Fermi energy £p of & pure atomic Fermi gas. We do not 
present the phase diagram for the molecular fraction rj m — N m /N tot , since 
it is very similar to that given by Williams et for an ideal gas. 



3. Adiabatic Phase Diagrams 

In this section, we discuss the adiabatic phase diagrams. We first briefly 
review the experimental procedures™' In these experiments, the ultracold 
two-component Fermi gas is initially prepared at a magnetic field detuned 
far from the resonance position, corresponding to a pure atomic gas. In 
this weakly interacting regime, the temperature of Fermi gas is measured 
using time-of-flight imaging. The magnetic field is then slowly lowered to 
the vicinity of the resonance to allow the atoms and molecules sufficient time 
to move and collide in the trap. This process may satisfy the condition that 
the gas be able to collisionally relax to equilibrium P We regard this process 
as an adiabatic quasistatic process. The system will follow a path of constant 
entropy in the phase diagram. 

In FigEfA), we show contours of constant entropy, which is the sum 
of the atom entropy and the molecule entropy, in the e res - T plane. We 
recognize those lines as the paths that are traversed as the resonance energy 
£ res is swept adiabatically. 

The relation between the initial temperature of the atomic gas and 
the final temperature of the molecular gas has been given for the ideal gas 
mixture of fermionic atoms and bosonic molecules by Williams et a/P This 
relation can be derived by connecting the initial entropy of the atomic gas 
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and the final entropy of the molecular gas. 

In FigEfB), we plot the final temperature in a deeply BEC region and 
a moderately BEC region as a function of the initial temperature. The solid 
line is obtained by a numerical calculation in a deeply BEC region where 
£res — — IOO^f- The dotted line is obtained by a numerical calculation in a 
moderately BEC region where £ res = — 4£p. The dashed line represents the 
high-temperature limit for an ideal gas, while the dot-dashed line represents 
the low-temperature limit for an ideal gas. The low-temperature approxi- 
mation agrees with the both of numerical solutions. The high temperature 
approximation does not agree with the numerical solution for a moderately 
BEC region, because atoms appear when the temperature is not sufficiently 
high and thus the concept of connecting the initial entropy of atoms and the 
final entropy of molecules breaks down. On the other hand, the high tem- 
perature approximation agrees with a numerical solution of a deeply BEC 
region except for the high initial temperature regime where atoms appear. 
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Fig. 2. Contours of constant entropy (A) and final temperature versus initial 
temperature (B). 

To compare our results with the phase diagrams of the experiments^! 
directly, we plot the fermionic condensate fraction r] c against the resonance 
energy e res and the initial temperature of a gas in FigEfA). The initial 
temperature T{ is found by equating the total entropy to the initial entropy 
Stot ( 5res ' ^) = Si(Ti). In FigEJ the transition temperature is T c ~ 0.537V 
in the BEC region where s res = — Aey- On the other hand, in FigEf A), the 
transition temperature in terms of the initial temerature is T^ c ~ 0.217V 
at £ res = — 4^p- This translated transition temperature T{ c in the BEC 
limit is comparable with the transition temperature observed in the resent 
experiments. ^1 

In FigOl^B), we plot the fermionic condensate fraction (black solid line) 
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and the condensate molecular fraction (red dotted line) with an initial tem- 
perature Ti/Tp — 0.08 as a function of the resonance energy. The blue 
dashed line represents the molecular fraction. The behavior of the conden- 
sate fraction shown in FigO^B) agrees with the experiments™! 
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Fig. 3. Adiabatic phase diagram (A) 
Ij = 0.08T F (B). 
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4. Molecular Conversion Efficiency 

Hodby et al w have investigated the molecular conversion efficiency 
They found that the molecular conversion efficiency does not reach 100% 
even for an adiabatic sweep, but saturates at a value that depends on the 
initial peak phase-space density. 

Assuming that the principle controlling the maximum molecule conver- 
sion efficiency proposed by Williams et al. A also applies when including the 
resonant interaction, we calculate the molecular conversion efficiency using 
Vm(^res — 0) as a function of initial temperature, as shown in Fig^I The red 
dotted line represents the result from our model, while the blue dashed line 
represents the ideal gas mixture of Fermi atoms and Bose molecules. The 
black dots represent the data of Hodby et a/.p^ and the black solid line rep- 
resents their simulation result. Black dotted lines represent the uncertainty 
of the pairing parameter in their simulation. Our model agrees well with the 
experimental data. 

The resonant interaction is more important for lower initial tempera- 
tures. At the initial temperature T{ — 0, the molecular conversion efficiency 
is 100% in the ideal gas mixture model. On the other hand, the experimen- 
tal data in the low temperature limit T{ — ► is suppressed from 100%. The 
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Fig. 4. The molecular conversion efficiencies. 

efficiency in our model at T{ — is also less than 100%, because the resonant 
interaction suppresses the molecular conversion. 

5. Conclusion 

The adiabatic phase diagram allows us to compare the theory with 
the experiments, since in the experientsP^l the bimodal condensate fraction 
have been measured against the initial temperature and the magnetic field. 
We obtain the transition temperature in terms of the initial temperature 
Ti C — 0.217V in the BEC regime, which is consistent with the experimental 
result PEI Finally, the molecular conversion efficiency is plotted as a function 
of the initial temperature. We find that the resonant interaction suppresses 
the complete conversion. 
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